Aging in a simple model of a structural glass 
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Abstract 



1: 

constraints in contact with a particle reservoir. Quench below the glass transition is represented 
by the jump of the chemical potential above a threshold. After a quench, the density approaches 
the critical density — where the diffusion coefficient of the particles vanishes — following a power 
. law in time. In this regime, the two-time self-correlation functions exhibit aging. The behavior 

of the model can be understood in terms of simple mean-field arguments. 

o 1 INTRODUCTION 

. Although glasses are considered the paradigmatic example of long-lived violation of thermodynam- 
i ical equihbrium, the detailed nature of this violation is still the subject of much debate [0, 0]. 
Phenomenologically, a glass is an undercooled liquid, quenched below a temperature such that its 
^ ' viscosity exceeds 10^^ Poise. In this situation, the glass can be considered "sohd" in the sense that its 
CN ■ molecules rattle within a "cage" formed by their neighbors and do not leave it, at least within exper- 



imental time scales. However, even in the glassy state, the system does not reach a time-independent 
Q\ i statistical state, but rather keeps evolving at a slower and slower pace as the time elapsed since 
its quench increases. This is the origin of the striking aging effects observed in glasses [§, 0. Aging 
corresponds to the property that, while one-time quantities like the average energy, volume, etc., 
appear to be invariant under time translations, two-time correlations and responses exhibit a non 
trivial dependence on both of their arguments [^]. Similar properties have also been observed in 
O ', spin-glasses and in fact appear quite frequently, provided one takes the infinite size limit before 
'p. ' taking the infinite time limit It now appears that aging can be described by a consistent and 
robust phenomenology, in which several apparently unrelated properties are brought together. On 
><| ■ the one hand, in purely relaxational systems, the lack of time-translation invariance (TTI) in the 
. correlation functions is associated with the violation of the fluctuation-dissipation theorem (FDT) 
" " ' [^. This violation can be related to entropy production and leads to the definition of an effective 
frequency-dependent temperature out of equilibrium |T^. 



Aging properties have been thoroughly investigated within spin-glass models [|TTI] which are closely 
related to the mode-coupling theories of structural glasses [|T], The spin-glass models proposed 
as a description of structural glasses lack, however, a transparent physical interpretation in terms of 
particles and involve a complex, random hamiltonian which is hard to justify as a description of a 
fluid. Their justification is rather a posteriori in the sense that the phenomena they exhibit recall 
the behavior of structural glasses. On the other hand a class of very simple kinetic models have been 
introduced to describe the slowing down of the dynamics []T2|, |13| . These models are defined by kinetic 



rules involving a selection of the possible configuration changes ("moves") and are therefore called 
models with constrained dynamics. The kinetic rules satisfy detailed balance and are compatible with 
a Boltzmann-Gibbs equilibrium distribution involving a hamiltonian, usually chosen to be a trivial 



one. The constraints are alone responsible for the slowing down of the dynamics because, near any 
allowed configuration, there are only few configurations which satisfy them. Since the hamiltonian is 
trivial, it is easy to prepare the system in an equilibrium state. By following its subsequent evolution, 
one can then collect data on the equilibrium correlation functions, which exhibit TTI by definition. 
In order to observe aging and FDT violation, it is necessary to prepare the sample in an out-of 
equilibrium state. One way to achieve this goal is to introduce a conjugate field which plays the role 
of the temperature, and to simulate the quench by a drastic change of this field. 



In the present contribution, we consider a lattice-gas model introduced by Kob and Andersen 14 



and generalize it to allow particle exchange with a reservoir |T3 . The kinetic constraints prevent the 



particle from moving when it has too many neighbors: the detailed balance condition is satisfied if 
the particle cannot move also if it would have too many neighbors after the move. As a consequence, 
the evolution of the system slows down and becomes sluggish when the particle density p increases: 



when p > pc ^ 0.88 the self-diffusion constant D vanishes ||l4 

We have allowed the system to exchange particles with a reservoir, introducing therefore the 
intensive variable p/k-oT (the ratio of the chemical potential to the absolute temperature), conjugate 
to the total number of particles. Since the hamiltonian is trivial, temperature does not play any 
significant role and we can set k-oT = 1 throughout. The equilibrium equation of state p = Pcq{p) 
can be then trivially calculated. There is therefore a critical value pc of p defined by Peq{pc) = Pc- 
This value plays the role of the glass temperature. In particular a quench is represented by a jump 
in p from below to above pc, corresponding to a sudden compression. 

Remarkably, this simple model exhibits a number of glassy properties so far obtained only in 
much more complicated models One can perform numerical experiments analogous to smooth 
cooling, quench, and hysteresis cycles: namely, one can let p increase smoothly, perform a cycle, 
or jump suddenly above pc- In the first case the density p first increases, then reaches a plateau 
value that depends on the compression speed: the slower the compression, the higher the value. The 
critical density pc could apparently be reached in the limit of zero speed. In the second case, the 
density appears to follow a hysteresis cycle whose area decreases as the compression speed decreases. 
And finally, in the case of a quench, the density appears to approach the critical one like a power 
law. In this case, aging effects in the self-correlation function B{t,t') = {\f{t) — r{t')\'^) are evident. 
We show here that the power-law relaxation of the density can be interpreted as a consequence of the 
fact that the diffusion coefficient vanishes at pc, and that this in turns explains the aging effects in 
B{t, t'). A simple mean-field-like theory is proposed, which could act as a guide for further numerical 
investigations. 

Section 2 contains the definition of the model and the results of its simulation. Sec. 3 contains a 
simple mean-field theory which accounts for the behavior of the model under quench. Sec. 4 contains 
a brief conclusion and hints for further investigation. Some numerical results reported in Sec. 2 had 



already been published in ref. [15 



2 NUMERICAL RESULTS 

We consider the kinetic lattice-gas model studied by Kob and Andersen (KA) ||14|. The system 



consists of particles in a cubic lattice of side L, {V = L^) with periodic boundary conditions. 
There can be at most one particle per site. Apart from this hard core constraint there are no other 
static interactions among the particles. At each time step a particle and one of its neighbouring sites 
are chosen at random. The particle moves if the three following conditions are all met: 

1. the neighbouring site is empty; 

2. the particle has less than m nearest neighbours; 



3. the particle will have less than m nearest neighbours after it has moved. 



The rule is symmetric in time, detailed balance is satisfied and the allowed configurations have the 
same weight in equilibrium. Significant results are obtained when the threshold m is set to 4 for a 
square lattice in three dimensions. With this simple definition one can straightforwardly proceed to 
study the dynamical properties of the model at equilibrium. 



Kob and Andersen have studied this model starting from an equilibrium configuration, which 
is obtained very simply by filling in turn each site of the lattice with a probability equal to the 
density p. In this way, TTI is satisfied by construction. One observes, in particular, that the 
diffusion coefficient of the particles vanishes above the critical density pc — 0.88 , and behaves like 
\p — Pc\'^ for p smaller, but close to pc- 

In order to study aging effects, we have generalized |jl5| the KA model by allowing the particles 



to appear and disappear on a single two dimensional layer (the "surface"). This event takes place 
with the following Montecarlo rule: a site on the surface is chosen at random; if it is empty we add 
a particle, otherwise we remove the particle with probability e~'^. We assume > 0, since however 
the interesting behavior takes place at high density (p > 0.5). Because of the periodic boundary 
conditions, the sample can also be conceived as a slab, whose upper and lower surfaces (placed, say, 
aX z = ±L/2) are in contact with the reservoir. For each value of p it is trivial to evaluate the 
equilibrium density 

p,^{p) = {l + e~^y\ (1) 

If Pcq(p) < Pc the system rapidly reaches the equilibrium state. There is therefore a critical value, 
/Xc, of p, defined by Peq(pc) = Pc- A quench corresponds to a jump of p from below to above pc. 
Therefore, p plays a role analogous to the inverse temperature in mode-coupling theories. We then 
observe (see fig. |l]) that after a quench p never exceeds Pc, but rather approaches it like a power law 
in time: 

p(r) oc r-^ (2) 

where r is the "effective time" after the quench (see later) and the exponent z is approximately equal 
to 0.3 .ITBIl Let us remark that zd) = 1 within the errors. 



The self-correlation function B{t,t') is defined by 

B{t,t') = {m-m\'), (3) 



where f{t) and f{t') are the positions of the same particles at times t and t' respectively. The average 
must be defined with some care, since the particles may leave or enter the system. We define it by 
averaging only over the particles which are present at both times. 

From now on we denote the smaller of the two times t, t' by (the "waiting time"). Aging 
corresponds to the fact that B(t,t^) does not depend only on t — t^, but rather appears to be a 
function of t/tw Similar results have been obtained for the p-spin spherical spin-glass model [|]. 
From the figure it is evident that for small values of TTI approximately holds, and aging sets in 
only for larger values of t^- In fact, we find here a phenomenon already encountered by Kob and 



Barrat [0 in Lennard- Jones systems: one has to consider an effective waiting time that takes into 
account the relaxation time of the system before the quench. Hence, one should define an effective 
waiting time Tw = tw + tq, where tq is the relaxation time characteristic of the equilibrium situation 
|T7| . Figure ^ shows that, if one plots B{t + t^,t^) vs. t/r^, the curves lie roughly on top of one 
another. The deviations are probably due to finite size effects, e.g., to particles that escape from the 
sample. Indeed, since the contribution of these particles to B{t + t„,tw) has not been included, the 
mean squared displacement has a tendency to be underestimated at long times. The plot exhibits 
the waiting-time dependence of the diffusion constant and is consistent with a simple argument 
suggesting that diffusion is logarithmic in the presence of aging ||1 
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Figure 1: Relaxation of the excess specific volume 5v — v — Voo, where v — 1/p, Voo — 1/Pc, after a 
quench to the subcritical value 1/ fi = 1/2.2. The density of initial configuration is 0.75. The straight 
line is Sv = 1.36 • t'^-^^^, with = 0.88 . 
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Figure 2: Mean square displacement vs. time after a quench to the value l//i = 1/2.2. The waiting 
times are = 10, 10^, 10^, 10"^, 10^ . The density of initial configuration is 0.75. Average over ten 
samples. 



3 MEAN-FIELD ARGUMENTS 



In this section we show how we can understand the results reported above by simple mean-filed-like 
arguments. Let us assume, as found by KA, that -D(p) ~ (pc — pY^ where, as found by KA, ~ 3.1 . 
We consider the sample as a slab, whose free surfaces sX z = ±L/2 are in contact with the reservoir 
and therefore rapidly reach a density equal to the critical one. Within the sample, the density p{z, t) 
satisfies the following non-linear diffusion equation: 

dp d 
dt dz 



(4) 



with the boundary conditions p{L/2) = p{—L/2) = p^, and where p{z,t) < p^ for —L/2 < z < L/2 
and all t. Changing variable y = p — pc "we obtain the equation 



One can look for solutions of the form y{z,t) = Y{t)f{z). One obtains 

Y{t) ~ r^/"^, (6) 

corresponding to a power-law relaxation of the density with an exponent z = 1/(1), and a differential 
equation for the density profile f{z): 

/(^) = §; [fH^)f{^)] , (7) 

with the boundary conditions 

/(±L/2) = 0. (8) 
It is easy to see that this equation allows for solutions of the form 

f{z) = hf{z/L), (9) 

where 

/oocL2/(^-<^), (1q) 
and where the universal profile f{x) satisfies 

/(x) = 0/>-i(/')^ + />r, (11) 

where one can choose the boundary conditions /(O) = f{0)/fo = 1, /'(O) = 0. The result of a 
numerical integration of this equation with the KA value of = 3.1 is shown in fig. ^. The solution 
vanishes for x = Xq 0.72 . The full solution is therefore f{z) = f (0) f {2xoZ / L) . Unfortunately our 
simulation data are still too noisy to allow for a meaningful comparison with this prediction. 

Let us now consider the effects of a varying diffusion constant on the self-correlation function 
B{t,t^). In a simple minded approach we would obtain 

B{t,t^) = f dt'D{t'). (12) 

If we assume D{t) oc we obtain 

5(t,t„)oc(ti-^-C^). (13) 
However, in our case, we have C, = zcf) = 1, and a simple integration leads to 

S(t,t^) oc (logt-logtw). (14) 

This is the functional obtained in ref. |T8| from the hypothesis that B{t,t^) depends only on t/t^ 
and from the "triangle relation" for B{t,t'): 

B{t,t') = B{t,s) + B{s,t'), ioTt'<s<t, (15) 

which stems out of the statistical independence of particle displacements over nonoverlapping time 
intervals. This functional form is borne out by the simulations, as can be seen by fig. H. 
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Figure 3: Universal density profile f{z)/f{0) vs. z/ {L/2). We have chosen the KA value of = 3.1 . 
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Figure 4: Mean squared displacement B{t + t^^t^) vs. scaled time t/r„ (where r„ = + tq) for 
= 10, 10^, 10^, 10^, 10^, and tq = 10^. Quench to the subcritical value = 1/2.2, starting from 
density 0.75. The full line corresponds to B{t + t^,t„) oc (log(t + t^,) — logtw)- System of size 20^, 
average over ten samples. 



4 DISCUSSION 



We find remarkable that this toy model — based as it is upon indefendable assumptions — exhibits a 
behavior so similar to those of much more complex models. Perhaps the most appealing aspect is 
the fact that aging appears as a consequence of the slow (power-law) approach to the critical density. 
This approach bears some similarity to self-organized criticality: after a sudden compression, the 
system endeavors to accommodate as many particles as it can: but, if the critical density is exceeded, 
particle diffusion stops and no more particles can get in. The result is an ever slower approach to the 
critical density, and it may be said that the system ages because it approaches criticality. Similar 
considerations in the context of the Bak-Sneppen model have been made by Boettcher and Paczuski 
p!9| . In our case they can be made explicite and quantitative, even within a mean-field-like approach. 
Several interesting informations concernig the breakdown of thermodynamical equilibrium can be 

— * 

obtained from response functions. In our case, one can consider applying a small constant force Fa 
to particle a between times tw and t. As a consequence, the particle suffers an extra displacement 
Sra(t,t„). The response function t^) is then defined as 



X{t,t^) = , (16) 

where i is one of the coordinates. The celebrated CK plot ^ of the response function vs. the corre- 
lation function should then appear as a broken line, where the break takes place in correspondence 
of the value of i?(t,t„) for which TTI no longer holds. If the "mean-field" arguments given above 
hold, the response of the aging system should be locked in at what it would exhibit at the critical 
density. The slope of the FDT violating line should then be independent of the quench value of /j,. It 
would be nice to be able to put this slope in relation with "universal" properties of the model, like 
6, but we have so far been unable to obtain this relation. 
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